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In two spatial dimensions quantum mechanical particles are not limited to 
being bosons or fermions as they are in three dimensions, but can be particles 
known as anyons yQ-S] which obey braiding statistics — meaning that the quan- 
tum mechanical state of a system is altered in a particular way when one anyon 
is moved (braided) around another, independent of the specific path chosen. 
Such anyons come in two major varieties — Abelian and non-Abelian — both of 
which were long ago predicted |2]-|5] to be realized in certain Fractional Quantum 
Hall (FQH) systems. Non-Abelian anyons have recently received considerable 
attention as a promising platform for topological quantum computation [i3j, |6]- 
[Sj. However, experimental demonstration of anyonic braiding properties has 
remained elusive and very controversial | [9| 110] . In this paper, we report on the 
study of Abelian and non-Abelian statistics through Fabry-Perot interferometry. 
Our detection of phase slips in quantum interference experiments demonstrates 
a powerful, new way of detecting braiding of anyons. We confirm the Abelian 
anyonic braiding statistics in the v = 7/3 FQH state through detection of the 
predicted statistical phase angle of 2n/3, consistent with a change of the anyonic 
particle number by one. The z/ = 5/2 FQH state [11] is theoretically believed to 
harbor non-Abelian anyons [2], [3] which are Majorana, meaning that each pair 
of quasiparticles contain a neutral fermion orbital which can be occupied or 
unoccupied and hence can act as a qubit. In this case our observed statistical 
phase slips agree with a theoretical modelfl^] where the Majoranas are strongly 
coupled to each other, and strongly coupled to the edge modes of the interfer- 
ometer. In particular, an observed phase slip of approximately vr is interpreted 
as a sudden fiip of a qubit, or entry of a neutral fermion into the interferometer. 
Our results provide compelling support for the existence of non-Abelian anyons. 
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The most promising systems for the reahzation of anyons are the FQH states in two- 
dimensional electron systems at low temperatures and high magnetic fields. Each FQH state 
is labeled by its filling fraction the dimensionless ratio of electron density to magnetic flux. 
The primary FQH states at v = 1/m, with m odd, are thought to be an imcompressible 
quantum fluid. Theory predicts [H |5] that the quasiparticle excitations of these states should 
be anyons with fractional charge q = e/m and fractional statistical phase 6 = n/m occurring 
whenever two such anyons are exchanged in a clockwise manner. While there is convincing 
experimental evidence in support of the fractional charge [T3HT6], experiments on the study 
of fractional statistics have been received with skepticism [9]. 

The FQH state at = 5/2 is considered to be the most promising system for the demon- 
stration of non-Abelian anyon braiding statistics. The quasiparticle excitations are expected 
to be anyons with charge g = e/4 and whose braiding statistics have both an Abelian part 
with 6 = tt/8 and a more complicated non-Abelian part which associates a Majorana with 
each quasiparticle [21 13]. 

Braiding statistics of anyonic quasiparticles for FQH systems may be detected through 
quantum interferometry experiments, in which test quasiparticles "braid" a group of local- 
ized quasiparticles at the center of the interferometer [31 [814T01 fTT l flQl - IM] . In a quantum Hall 
Fabry-Perot interferometer, interference trajectories are realized through coherent transport 
along edge channels and tunneling across a pair of split-gated constrictions, which act as 
beam splitters (See Fig. lb). In the Abelian case, assuming that tunneling at the first and 
the second constrictions occur with amplitude ti and t2, the conductance, G, across the 
interferometer is 

Goc |tip + |t2p + 2Re{ttt2e^'^} (1) 

The resulting interference phase, 0, is a sum of the Aharonov-Bohm phase due to the 
magnetic flux enclosed by the edge state trajectory and the statistical phase 29 (a loop 
equals two exchanges, hence 26) for each quasiparticle encircled by the edge. The general 
scheme of this measurement is to observe discrete phase-slip events of 26 associated with 
the sudden entry or exit of a quasiparticle [T2l [25] from the interferometer. 

The Fabry-Perot interferometer was fabricated from a high mobility, symmetrically doped 
GaAs/AlGaAs quantum well. The layout of the interferometer is shown in Fig. la. The 
diameter of the interferometer is 1.2/im, and the width of the two constrictions is ~400 
nm. The sample was mounted on a dilution refrigerator capable of reaching below 10 mK. 
[Magnetotransport through the interferometer is shown in the supplementary material.] The 
diagonal resistance, R^, is measured as the voltage bias on the plunger gate (lower middle 
gate of Fig. la) is swept at a steady rate. The bias changes the encircled area resulting in 
interference oscillations. Experimental studies of interference in the integer and fractional 
quantum Hall regimes have been reported earlier [HI UHl I^TH2^ . The effect of Coulomb 
charging in the interferometer may be either strong or weak, in which case we say the device is 
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in either the Coulomb-dominated or Aharonov-Bohm regime respectivelyj21l |23| IMf [26| |27] . 
The phases observed (for example the gate period) can be modified by such Coulomb effects 
if they are present. However, it has also been discussed [12] that under certain reasonable 
conditions, the phase slip when a quasiparticle enters the interferometer could be immune 
to these modifications even in the Coulomb-dominated regime. 

Fig. 1 illustrates sweeps of the plunger gate voltage Vp for the z/ = 2, 7/3 and 5/2 
FQH states. In the case of the v = 2 state (Fig. Ic), Rd oscillates periodically with Vp. 
For the 7/3 state (Fig. Id), Rp) shows a periodic behavior for Vp < —27 mV. For —27 
mV <Vp < — 18mV, the behavior shows a telegraph noise behavior between two curves 
(black and blue) that are out of phase by exactly A0 = 27r/3. This noise refiects the phase 
associated with an anyon randomly entering and exiting the interferometer. Another 27r/3 
phase slip can be identified near Vp = — 15mV [See Supplementary Material for numerical 
analysis]. Observation of these phase slips demonstrates the predicted braiding statistics of 
the e/3 anyons for the 7/3 FQH state. 

For the 5/2 state, the above equation for the interference may be modified to account for 
the non-Abelian properties of the quasiparticle excitations [3] [8| [T81 - I20] . In this case, if there 
is an even number of quasiparticles in the interferometer, interference should be observed, 
but the phase of the interference may be shifted by vr depending on the parity of the number 
of neutral fermion orbitals that are occupied inside the interferometer (i.e., the settings of 
the qubits). A fiip of a qubit can suddenly occur (equivalent to a neutral fermion entering 
the device) resulting in a vr phase slip in the interference |12j. 

If there are an odd number of e/4 quasiparticles inside the interferometer, no interference 
is predicted so long as the quasiparticles are isolated from the edge[T9l |20]. However, for a 
small enough device, quasiparticles will always be strongly coupled to the edge. We believe 
that our device is in this limit. In this case [121 ESI EH], instead, it is predicted that the 
non-Abelian part of the quasiparticle should be absorbed into the edge; interference does 
not vanish, and one observes only the phase 26 = tt/A from the Abelian part. Thus the full 
prediction [12] for this device is that one should observe vr phase slips associated with fiips 
of the qubit, and 7r/4 slips associated with quasiparticle addition (both can occur together 
resulting in a slip of 57r/4). 

Fig. le illustrates a sweep of the plunger gate voltage and the analysis of the phase 
slips for the 5/2 FQH state. Somewhat similar to the 7/3 data, Rpi shows an approximate 
periodic behavior with a clear phase slip occurring around -35 mV. The measured phase slip 
angle is very close to A0 = 1.257r in agreement with theoretical prediction. 

Fig. 2a shows the detailed interference profile of the z/ = 5/2 FQHE state for many 
plunger gate sweeps at 21 successive values of magnetic fields. (The general structure of the 
plot indicates that this device is Coulomb-dominated). A striking feature is that R^ does 
not change monotonically with gate voltage and magnetic field. There are instances where 



3 



a minimum becomes a maximum and vice versa (i.e., phase slips of approximately vr). 

In Figs. 2b-2f show representative plunger gate voltage sweeps from Fig. 2a. Fit of the 
Kd sweeps show that an overwhelming majority of the sweeps exhibit clear phase slips of 
A0 vr to 57r/4 in the range of gate voltages between -32 and -37 mV. 

In order to make a unbiased analysis of these phase slip angles, we have developed a 
computerized algorithm to best determine the phase slips present in Fig. 2a. (See Supple- 
mentary Material.) At each point in the {Vp,B) plane, the gate sweep is fitted with two 
sine wave segments, one for V" > Vp and one for V < Vp, and the relative phase slip between 
these two curves is plotted in Fig 3a. The red background represents the phase shps near 
and 27r (meaning no phase slip). The nontrivial regions of this plot are then histogrammed 
in Fig 3b. Two prominent persistent features are observed: most clearly, strong phase slips 
ranging from tt to 57r/4 are observed in the range of gate voltages from -42 to -31 mV (blue 
to light blue). Phase slips near .357r (yellow-green) are observed between -27 to -17 mV. 
Much more weakly, both the yellow-green and blue regions show slips near O.Ttt. 

Given that our system is Coulomb-dominated, the theoretically expected phase slip of 
IT /A does not have to hold precisely [12], and it is not unreasonable to conclude that the 
measured 0.357r is a reflection of a 7r/4 slip with some modification. The weak O.Ttc peak 
could be a pair of 0.357r slips in a row. It is also possible that the slips between vr and 57r/4 
are one or the other value and these are either modified slightly by Coulomb effects, or are 
simply measured imperfectly due to other noise in our data. Whether or not we can pin down 
the expected ir/A and 57r/4, it is quite clear that slips very near vr are extremely dominant. 
While it has been pointed out^30j that interferometry experiments may not uniquely identify 
a particular quantum Hall state, such slips are clear evidence of neutral fermions, and are a 
crucial step towards demonstrating the reality of non-Abelian anyons. 

In summary, we have detected the braiding statistics of anyons. For the u = 7/3 case 
our data clearly demonstrates the predicted braiding statistics. For u = 5/2 our data is in 
agreement with current theory, clearly observing a neutral fermion. Finally, we note that 
it is quite difficult to obtain data of this type since too much telegraph noise can easily 
make interpretation impossible, and too little removes the desired slips. Nonetheless, with 
persistence this is a powerful technique for probing the anyonic properties of FQH states. 
We anticipate that the fidelity of anyonic detection will improve with time, possibly leading 
to the realization of a topological quantum information processor. 
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and W.K. acknowledge support by Microsoft Project Q and NSF MRSEC Program through 
the University of Chicago Materials Center (DMR-0820054). S.H.S. was supported by EP- 
SRC grants EP/I032487/1 and EP/I031014/1. L.N.P, K.W.W., and K.W.B acknowledge 
partial support from the Gordon and Betty Moore Foundation and the NSF MRSEC Pro- 
gram through the Princeton Center for Complex Materials (DMR-0819860). 
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Fabry-Perot interferometer and interference measurements. 

a, A SEM view of the Fabry-Perot interferometer (FPI). The low temperature mobihty 
of the unprocessed sample was 2.83 x 10''cm^/Vs with an electron density of 3.3 x 10^^ 
cm~^. The interferometer was initially defined using e-beam lithography, dry etched to the 
depth of the two-dimensional electron system, and metallized with TiAu prior to lift-off. 
The diameter of the interferometer is 1.2/im, and the width of the two constrictions is 
~400 nm. b. Expected edge state trajectory (red) in the quantum Hall regime for the 
FPI. Two narrow constrictions serve as beam splitters and interfere particles tunneling at 
the two constrictions. The observed interference phase is a sum of the Aharonov-Bohm 
phase due to the enclosed magnetic flux $ and the statistical phase due to the exchange, 
c. Diagonal resistance, Rd, as a function of the plunger gate voltage at the z/ = 2 quantum 
Hall state under B = 7.5614 tesla. Bias of the side gate from a changes the area enclosed 
by interfering particles and leads to a periodic oscillation, d. Measurement of Kn at the 
u = 7/3 fractional quantum Hall state under B = 5.7139 tesla. In conjunction with phase 
slips of A(j) = 2tt/3, telegraph noise is observed in some regions of gate voltage. Telegraph 
noise is most pronounced between -20 and 25mV, but telegraph noise behavior is seen above 
-20mV as well, presumably from discrete switching of the number of quasiparticles in the 
interferometer. The black curve is the best fit of the large bias portion of the data below 
-20mV. The blue curve is the same curve as the black curve but phase-shifted by 27r/3. The 
light blue is further phase shifted by 47r/3 relative to the black curve. With the phase slips 
occurring, the data can be parameterized from the black to blue to the light blue curves 
with the successive phase slips that are very close to 27r/3. e. Measurement of R/5 at the 
u = 5/2 fractional quantum Hall state under B = 5.5275 tesla. The blue curve is the same 
curve as the black curve but phase-shifted by 57r/4. Phase slip events occur between -35.5 
and -34 mV. 
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Figure 2: Diagonal resistance and phase slips at the u = 5/2 fractional quantum Hall state, 
a, The diagonal resistance, Rd, with a smooth background subtracted as a function of 
magnetic field, B, and plunger gate voltage, Vp. Rjj exhibits a series of oscillations whose 
minima and maxima tend to evolve with a positive dVp/dB curvature over the measured 
ranges of gate voltage and magnetic field. The dVp/dB curvature places the interference 
in the Coulomb dominated regime of interferometry where Coulomb interaction plays an 
important role, b-d Representative phase slips of vr and 571/4 between -37 and -30mV. The 
phase slip of ~ vr leads to a reversal of the parity of the Rj^. b: 5.5330 telsa, c: 5.5340 tesla, 
d: 5.5355 tesla. e. An example of sweep under B = 5.5315 tesla where there is a multiple 
phase slip events with A0 ^ ±57r/4 at -34mV and ~ ±7r/4 phase slips between -19 and 
-17mV. f. Plunger voltage sweep under B = 5.5305 tesla with multiple phase slip events with 
A(f) ^ ±71 between -35 and -33.5mV and ~ ±7r/4 phase slips between -23 and -18 mV. 
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Figure 3: Phase slips for the v = 'b/2 fractional quantum Hall state. 

a, Plot of phase slip angles as a function of plunger gate bias voltage and magnetic field 
from the diagonal resistance data from Fig. 2a. A computerized algorithm was used to 
determine the best phase slip in an unbiased way. b, Histogram of non-trivial regions of the 
phase slip angles from Fig. 3a. 
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Supplementary Fig. 1 compares the Hall resistance, Kxy, of the bulk, unprocessed part of 
the sample and the diagonal resistance, R^i, through the Fabry-Perot interferometer (FPI) 
in the second Landau level. These data were taken after an illumination by a red LED. 
Kxy in the bulk part of the sample shows fractional quantum Hall (FQH) states at filling 
factors u = 5/2,7/3,8/3, and 14/5. In addition, a pair of reenetrant insulating states to 
u = 2 and 3 are respectively observed on high and low sides of the u = 5/2 plateau. The 
FQH states within the FPI is noticeably weaker compared to that found in the bulk. The 
Kd through the interferometer exhibits a noticeable decrease in the electron density and 
degradation of various quantum Hall states. Only u = 5/2 and 7/3 states are visible and 
reentrant insulating states to z/ = 3 are suppressed. Such changes may be expected due to 
the small size of the interferometer and the side-depletion coming from the confinement. 
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REPRODUCIBILITY OF ANYONIC PHASE SLIPS 



An important question regarding our results of anyonic braiding and phase slips is how 
reproducible the data is. It tiu^ns out that observation of phase slip for both z/ = 7/3 and 5/2 
fractional quantum Hall states is rather difficult because it is not easy to realize the condition 
within the FPI under which phase slips can be detected via DC transport measurements. 
In order to realize the = 5/2 FQH state inside the FPI, we must illuminate the FPI using 
a LED during the cooldown. Too much or too little illumination can make it hard to detect 
any interference. Too much light activates low-frequency charge noise, which can overwhelm 
the interference signals. Too little light leads to low mobility device with no 5/2 state being 
realized and even no discernible interference signature. 

In our experience in working with many interferometers, whether or not a clear quantum 
interference can be detected depends critically on the electronic state that is realized within 
the two narrow constrictions that define the openings of the FPI. It is not always predictable 
whether or not the constrictions will behave like good beam-splitters. A large density differ- 
ence between two constrictions seems to be universally bad in suppressing the interference 
signal. With a reliably good interferometer - where good is defined as being able to detect 
interference signal and the FQH states in the second Landau level - we aim to optimize the 
illumination and cooldown procedure to realize uniformity between two constrictions and a 
low disorder FPI that can sustain robust FQH states. 

Another variability in the experiment is the timescale of the phase slip events, which is 
expected to be determined by the local potential. If the phase slips and the corresponding 
telegraph noise occurs at too fast of a time scale, then it becomes difficult to detect them 
via DC transport. If the time scale is too long, then it becomes difficult to collect sufficient 
number of events within a practical time scale to be feasible. 

With all these factors in play, it is not possible to detect clear phase slip events in 
interference measurements for every cooldown. However, we have succeeded in detecting 
phase slips in the z/ = 5/2 state on at least three separate cooldowns. In the two out of 
three cooldowns, we also detected phase slips in the v — 7 /2> state. (In the third cooldown, 
an equipment failure prevented the search for phase slip events for the v — 1 j?) state.) In 
Supplementary Figs. 2 and 3 we provide data of phase slips for the v — 1 /?> and z/ = 5/2 
respectively from the cooldowns #1 and #2. The = 7/3 phase slip data shown in Fig. 1 
of the main paper was taken during the cooldown #1. The v — b/2 phase slip data in the 
main paper came from the third cooldown. The phase slip angles and the best fit to the 



11 



data were verified by the computerized algorithm described in Section III. 
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Supplementary Figure 2: Phase slips for the v = 7/3 (upper) and v = 5/2 (lower) fractional 
quantum Hall states from cooldown ^1. For the v = 7/3 state, A(/> = 27r/3 was observed between 
-26 and -22mV. For the v = f)/2 state, A(/) ~ vr was observed at -31mV. The = 7/3 data in the 
Figure 1 of the paper came from this cooldown. 
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Supplementary Figure 3: Phase slips for the v = 7/3 (upper) and v = 5/2 (lower) fractional 
quantum Hall states from cooldown #2. For the v = 7/3 state, Ac/; = 27r/3 occurs at -12.5mV. 
For the v = 5/2 state, A^ vr phase slip occurs multiple times between -12.5 and -9mV. 
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PHASE SLIP FITTING ALGORITHMS 



Slips are identified separately for each sweep of the gate voltage Vp. For each point Vp 
the program maximizes a fit quality of the measured resistance Rp, to two model sine-waves. 
We propose sine wave TZ]^ for Vp < Vp and sine wave TZp) for Vp > Vp. The form of these 
proposed model functions are 

7^+(yp) ^a + bVp + ccos(fVp - 0+) 
7^-(\/p) ^a + bVp + ccos{fVp - (/)_) 

where a, b, c, f, 0+, 0_ are all free parameters (note that a linear drift term b is also included). 

The fit is performed over the range Vp — 6V < Vp < Vp + 6V for some fixed range SV. 
The fit quality that is maximized is the following function 

Q= / dVpF[RD{Vp)] 

where 

f[Rd{Vp)] = e{Vp - v^)g[Rd{Vp) - n+{Vp)] + e{v^ - Vp)g[Rd - {Vp)n^{Vp)] 

Here the G is a nearness function for the fit 

^[""^ " 1 + Jx/SW 

with SR a fixed parameter. The nearness function is maximized if the model function is 
exactly equal to the measured value (i.e., for x = 0). Essentially this is least squares fit, but 
the penalty is softened at larger deviation from the model > 6R) so as not to heavily 
penalize outlier points which are very far from the model. Above O plays the role of the 
Heaviside step function, so that one tries to fit the measurement to the model function TZ'^ 
for Vp > Vp and to TZ]^ for Vp < Vp. For better numerical performance of our algorithm 
and smoother results is smeared somewhat and instead takes the form 



Q{x) 



1 x>0 
e"!"" X < 



with some fixed smearing distance w. Fixing the three w, 51^, and 5/2, the remaining 
parameters a,b,c, f,4>+,4>- are optimized numerically to maximize the quality of the fit 
Q. Since this is a nonhnear maximization, occasionally the algorithm (Nelder-Mead) will 
numerically converge to a local rather than global maximum. To combat this problem, the 
algorithm is run ten times with different seed values and the best result is chosen (typically 
only a few runs will not converge to almost the same result). To remain unbiased, the phases 
0+ and 0_ for the seed values are always chosen randomly. 
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Once the global maximum quality fit is found, the measured phase slip at the point Vp 
is then given by 0+ — </>_. Results do not differ much for a wide range of values of the three 
chosen parameters w, 6V and 6R. 
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Supplementary Figure 4: Best fit phase slip as a function of gate voltage &t v = 7/3. For 
most values of the gate voltage the measured slip is very close to or 27r implying no phase slip 
is detected. However, over two large regions phase slips are detected very close to 47r/3 (marks in 
bold). The blue points are the positions in gate voltage where fits are explicitly displayed in Fig. 5 
below. 

In Supplementary Fig. 4 we show the fit value of the phase slip 0+ — 0_ for each value 
of Vp for the data taken at u = 7/3. (In the fitting program, 6R = .OlA;^, w = .5mV and 
6V = 4 mV). At most values of gate voltage the measured phase slip is very close to or 
27r indicating no phase slip. However, over two regions of gate voltage the best fit clusters 
around the value of An/ 3 (points are bold in the figure). The mean phase slip given by the 
cluster of (bold) points around Vp = —21mV is l.SOvr with a standard deviation of 0.047r. 
The mean phase slip of cluster of points around -15.5mV is 1.27 vr with a standard deviation 
of 0.047r. (If the slight outlier point on the right is removed from this set, then the mean 
goes down to 1.26 tt with standard deviation 0.027r. 

To show examples of the quality of these fits, we consider the two values of Vp marked in 
blue in Supplementary Fig. 4 and display the fits explicitly in Supplementary Fig. 5. The 
fits for u = 5/2 are made using a similar technique. (In this case we use 6R = .02k fi, 
w = ImV and 6V = 7 mV). Plots such as that in Supplementary Fig. 4, are made for each 
value of magnetic field and then these plots are re-assembled into Fig 3a of the main text. 
In Fig 3a one might see three notable features. (1) The prominent blue strip between -42 
mV and -30.7mV (2) The light magenta region between -30.7mV and -27mV, and finally 
(3) the yellow-green region between -27mV and -17mV — restricted to B values between 
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approximately i?=5.5275T and 5.5326T. We then "cut" these three rectangular regions out 
of the 2d plot 3a, and histogram all of the phase shifts in these regions, which are then 
shown in Fig 3b. For the blue strip, and to a lesser extent, the yellow region, peaks are very 
evident. For the magenta strip there is no clear peak in the data. 
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Supplementary Figure 5: Example fits. (Top) fit for = -21.2mV (Bottom) fit for = 
— 15.8mV. These two sample gate voltages are marked in blue in the above figure. In both top 
and bottom TZ'^ is marked as the thick part of the red curve, and Tl~^ is the thick part of the pink 
curve. The point where the curves go from thick to thin is exactly Vp, where the phase slip is 
being measured. 
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